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ABSTRACT 



Aims. We present the results of several multi-station Very Long Baseline Interferometry (VLBI) experiments conducted with the ESA 

spacecraft Venus Express as a target. To determine the true capabilities of VLBI tracking for future planetary missions in the solar 

system, it is necessary to demonstrate the accuracy of the method for existing operational spacecraft. 

Methods. We describe the software pipeline for the processing of phase referencing near-field VLBI observations and present results 

of the ESA Venus Express spacecraft observing campaign conducted in 2010-201 1. 

Results. We show that a highly accurate determination of spacecraft state-vectors is achievable with our method. The consistency of 

the positions indicates that an internal rms accuracy of 0.1 mas has been achieved. However, systematic eff'ects produce off'sets up to 

1 mas, but can be reduced by better modelling of the troposphere and ionosphere and closer target-calibrator configurations. 

Key words, astrometry - techniques: interferometric - instrumentation: interferometers - instrumentation: miscellaneous 



1. Introduction 

The VLBI spacecraft tracking technique has been success- 
fully applied to a number of deep-space missions. The tech- 
nique imposes minimal requirements on the mission's on-board 
instrumentation and can operate with almost any radio sig- 
nal from a spacecraft, provided it is sufficiently powerful and 
phase- stable. Among VLBI-tracked missions one can mention 
the VEGA Venus atmosphere balloons JPreston et al. 1986 ), 
the Huygens Probe during its descent to the surface of 
Saturn's moon Titan (Lebreton et al. 20031) , VLBI tracking with 
the European VLBI Network (EVN) antennas of the con- 
trolled impac^_o£^SA^s_Smart-l probe on the surface of the 
Moon ( [Pogrebenko et al. 20Q6| ), VLBI tracking of NASA's Mars 
Exploration Rov er B spacecraft during its final cruise phase 
( |Lanyi et al. 2007 J, VLBI tr acking of the Cassini spacecraft 
at Saturn JJoneset al. 20111) , and the latest results achieved 
with ESA's Venus Express (VEX) VLBI spacecraft obser- 
vations and the ESA M ars Express (MEX) Phobos-flyby 
dMolera Calves et al. 201 Ob monitored by the EVN radio tele- 
scopes. 

The Planetary Radio Interferometry and Doppler 
Experiment (PRIDE) is an initiative by the Joint Institute 
for VLBI in Europe (JIVE). Designed as a multi-purpose, 
multi-disciplinary enhancement of planetary missions science 
return, PRIDE is able to provide ultra-precise estimates of 
spacecraft state vectors based on the phase-referenced VLBI 
tracking and radial Doppler measurements. These can be used 
for a variety of scientific applications including planetary sci- 
ence (measurements of tidal deformations of planetary moons, 
atmosphere dynamics and climatology as well as seismology, 



tectonics, internal structure and composition of planetary 
bodies), ultra-precise celestial mechanics of planetary systems, 
gravimetry and fundamental physics (e.g. tests of General 
Relativity or study of the an omalous accelerations of deep- space 
probes ( [Turyshev & Toth 2Q10| )). Spacecraft Doppler tracking is 
known to be the only possible wa y of detecting low-frequency 
(10"^ - 1 Hz) gravitational waves ( [Kopeikin & Schaefer 1999 ). 
Studies of interplanetary scintillations (Molera Calves et al., 
in prep.) are among the potentially rewarding by-products of 
PRIDE. PRIDE can also provide diagnostics of deep-space 
missions ("health check") and direct-to-Earth delivery of a 
limited amount of critical data (e.g. Friedman et al. (l2QQ8b ). 



PRIDE can be conducted with literally any radio- 
emitting spacecraft. Current PRIDE 'customers' include 
ESA's VEX and MEX spacecraft and th e Russian Spac e 
Agency's Space VLBI mission RadioAstron JKardashev 2Q09b . 
PRIDE has been consid ered in a number of design stud- 
i es, such as E VE JWilsonet aL 2011b, JUICE/Laplace 
(Blanc et al. 2009 , [Christophe et al. 2QQ9|), Ta ndEM 

( Coustenis et al. 2QQ9J , Kronos ( [Marty et al. 2QQ9| ), and 



several others. 



In this paper, we describe the principles and pipeline of the 
data acquisition, processing and analysis within the scope of 
PRIDE and present the results of the first large-scale observa- 
tional EVN campaign conducted in 2009-201 1 with ESA's VEX 
spacecraft. 
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2. Observations and data processing pipeline 

In phase-referencing VLBI, the relative position of a specific 
source is determined by observing a known source nearby (i.e. at 
a short angular distance on the celestial sphere), called a phase- 
reference calibrator, and applying calibration values of its fringe 
phase to the fringe phases of a target source. The most accurate 
measurements can be achieved when both the target and cali- 
brator sources are within the same primary beam of a VLBI 
antenna or no more than a few degrees apart ( Ros et al. 1999i) . 
VLBI phase referencing for spacecraft tracking involves alter- 
nated observations of the radio signal from the spacecraft and 
a nearby calibrator, most commonly a quasar. The best switch- 
ing time between the calibrator and target sources for out-of- 
beam phase-referencing ranges from tens of seconds to several 
minutes. In the out-of-beam regime, the ideal scan duration on 
each source depends primarily on the efficiency of the radio tele- 
scopes, their system temperature, and the slewing speed. In the 
observations described below, the scan duration was set to 2-5 
minutes. 

Spacecraft VLBI-tracking data processing and analy- 
sis tools are based on two major software engines: the 
ultra-high spectral resolution and phase detection soft- 
ware package SWSpec/SCTrack/Phase-Lock-Loop, developed 
at the Metsahovi Radio Observatory in collaboration with 
JIV^ and the EVN Software Correlator at JIVE SFXC 



Raw observational data 



(Keimpema et al. 201 1 ). 

Successful detection, processing, and analysis of space- 
craft VLBI-tracking data requires usage of a near-field delay 
model and computation of near-field analogues of baseline uv- 
projections. To address these tasks, a software package VINT 
(VLBI in the near-filed toolkit) was developed at JIVE and is 
described in detail below. 



2.1. Observations 

Preparation for observations, scheduling. In standard VLBI 
observations, coordinates of a source are considered to be con- 
stant during the observing run for all participating telescopes. 
Usually, the source position in the geocentric equatorial inertial 
coordinate system J2000.0 is used. For near-field VLBI obser- 
vations of spacecraft, the target can move rapidly across the pri- 
mary beam of a telescope. Therefore, the position of a spacecraft 
is calculated at diff'erent epochs and is then used for pointing the 
antenna at those particular epochs. In the "ultra" near-field case, 
i.e. when the target is within (or very close to) the synthesised 
aperture of a VLBI arra}Q so-called "cross-eyed" scheduling 
should be used. This means that the spacecraft J2000.0 positions 
are diff'erent in the schedules for diff'erent VLBI stations. 

Raw data. The down-converted signals observed by ra- 
dio telescopes are recorded with the standard VLBI record- 
ing systems Mark5 A/B (Whitney 1999) , or the PC-EVN 
( [Mujunen & Ritakari 2QQ4| ). Both systems were designed exclu- 
sively to record formatted data of astronomic and geodetic VLBI 
observations, allowing high data-rate ranging up to 2 Gbps. 

Fig. [T] represents the data processing pipeline that is de- 
scribed in detail below. 



and Pogrebenko, S.V. 

Spectrometer and 

Release, GNU GPL, 



^ Wagner, J., Molera Calves, G. 
2009-2012, Metsahovi Software 

Spacecraft Tracking tools, Software 
http://www.metsahovi.fi/en/vlbi/spec/index 

^ this is true for VLBI observations of, e.g., GNSS satellites or the 
Russian Space Agency's orbital radio telescope RadioAstron 
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Residual group delay 
and phase 



Residual phases of the 
carrier and ranging tones 



Delay/phase corrections 



Group and phase delay of the S/C signal 
with resolved 2k ambiguity 



Reconstruction of the apparent 
state vectors of the S/C 



Fig. 1. Data processing pipeline of phase-referencing VLBI ob- 
servations of spacecraft. 



2.2. Broadband correlation of reference-source signal and 
S/C data band 

Initial processing. Once at JIVE, the raw data are 
correlated with the JIVE Software Correlator SFXC 



(Keimpe ma et al. 2011| ), which has been developed at JIVE and 
is currently used as a production correlator for some EVN obser- 
vations. The SFXC is based on the original design developed for 
VLBI tracking of the Huygens Probe (P ogrebenko et al. 2008| ) 
and is capable of supporting both the far-field and near-field 
delay models and processing data in various "flavours" of the 
Mk5 format. 

When correlated with the SFXC, the fringe findeij^data are 
used for the clock search and bandpass correction. Also, a two- 
way Doppler phase correction (obtained with the SCTrack soft- 
ware, see next subsection) is applied to the spacecraft data to 
avoid frequency smearing. 

The correlation results are saved in a measurement set for- 
mat ( Ke mball & Wierin ga 2000). After this, the amplitudes and 
phases of the cross-correlation function are extracted for further 
processing. 

Group delay extraction. The complex cross-correlation 
functions (CCF) corrected for clock/clock-rate ofl'sets are 
smoothed in the frequency domain. After this, the phase is ex- 
tracted for each averaged spectrum and fitted by a linear func- 
tion. A group delay for the centre of each averaging interval is 
derived as the slope of the function Tgr = dcpf^/df, where / de- 
notes frequency. 

2.3. Narrow-band correlation of S/C signal 

SWSpec, SCTrack and phase-lock loop. Narrow band cor- 
relation of the spacecraft signal is performed with the 
SWSpec/SCTrack software package described in Molera Calves 
etal. (2011). 

The high-resolution spectrometer software (SWspec) allows 
the initial detection of the carrier tone of the spacecraft and tem- 
poral evolution of its frequency over the entire scan. In each 
SWspec pass, one of the raw data channels is extracted from 
the input data file. The software performs an accurate window- 
overlapped discrete Fourier transform (DFT) and time integra- 
tion of the spectrum. The generated time-integrated spectra are 
processed to determine the moving phase of the S/C tone fre- 
quencies along the scan. The series of frequency detections of 



a strong compact source observed to find interferometric fringes 
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the carrier line are fitted to a sixth-order phase-stopping polyno- 
mial. 

The core of the spacecraft tone tracking software (SCtracker) 
performs a phase polynomial correction to compensate for the 
Doppler shift, single and/or multi-tone tracking, signal filtering 
for each tone and, finally, phase detection of the carrier line with 
respect to the local H-maser clock. 

The final post-analysis of the narrow-bands is performed 
with the phase-lock loop (PLL) software. The PLL runs high- 
precision iterations of the phase polynomial fit on the filtered 
complex narrow-band signals. The implementation of the PLL 
allows filtering down the tones as often as desired to achieve 
mHz accuracy around the spacecraft carrier line or other tones in 
the band. The final residual phase in a specific stopped band is 
determined with respect to the set of subsequent frequency/phase 
polynomials initially applied. 

Reduction of phases to a common phase centre. 
Topocentric phase detections for each station are reduced to the 
common phase centre (usually the geocentre) by applying the 
geocentric signal delay, Tgc, for the near-field case (detailed de- 
scription of the model used is given in the next subsection): 



0gc(O = 0tc(^ - Tgc). 



(1) 



These geocentric phases are used as Doppler phase corrections 
for the broad-band correlation with SFXC. Finally, for each 



baseline i-j, diff'erential phases (ptj = i 



V,gc 



are formed. 



2.4. Resolution of the In-amblgulty 

A straightforward correlation of the wide-band spacecraft sig- 
nal allows deriving the group delay Tgr = dcp/df, while the 
narrow-band correlation of the S/C carrier tone allows deriving 
the phase delay Tph = 0//o, where /o is the observational fre- 
quency. These are two independent measurables, which should 
give the same phase solution. The rms error of the group delay 
is (T(Tgr) = l/(AF X S/N), with AF being the bandwidth of, e.g. 
8 or 16 MHz for the calibrator and about 1 MHz for the S/C but 
with an uneven distribution of S/N (signal-to-noise ratio) over 
the signal band. The rms error of the phas e del ay for a narrow- 
band signal is given by (T(Tph) = l/(/o x ^/S/N), where the S/N 
level for the spectral line is very high, so that even the ^/S/N 
may be larger than the S/N for the broad band case. In addition 
to this, /o » AF, so that the rms error of group delay is signifi- 
cantly larger than that of the phase delay. This in turn results in 
a diff'erent scatter of the phase solution resulting from group and 
phase delays. 

The group delay, although not that accurate, is a "true delay", 
while the phase delay is derived from the measured phase and 
may be subject to a In ambiguity. Since phases are 27i uncertain, 
the phase delay may have a bias of several cycles of ~ 120 ps for 
observations at X-band. Checking the phase delay versus group 
delay consistency by computing the uncertainty of the fit of the 
phase delay to the group delay and shifting the former by ^- 2;: 
(where n is an integer number) to obtain the best fit to the latter 
allows one to resolve the 27r ambiguity. 



2.5. Delay model 

The VLBI delay model is formulated in the bary centric celestial 
reference system (BCRS). This means that the delay is computed 
in the time scale of the BCRS and is subsequently transformed 
to the time-scale used for timing the measurements on Earth. 



Thus several time-scale and station coordinate transformations 
are required. 

The models described below are coded in Matlab and aggre- 
gated within the VINT software package. 

Geocentric terrestrial to celestial frame transformation. 
As a first step, a transformation of station coordinates from the 
terrestrial reference system (realised by the international terres- 
trial reference frame, ITRF) to the geocentric celestial reference 
system (GCRS) is conducted according to the latest recommen- 
dations of t he international Earth rotation service lERS 2010 
conventions JPetit & Luzum 20101) within the framework of the 
lAU 2000/2006 resolutions. 

Station positions and velocities in the ITRF at the epoch to 
(e.g. 1997.0 or 2000.0, depending on the realisation used), are 
reduced to the epoch of observation, taking the tectonic plate 
motion into account. A simple linear model is used: 



riTRF = To + To- (t - to). 



(2) 



At the next step, the station position, velocity, and accelera- 
tion in the GCRS are computed as follows: 

TgCRS = ^(0' TlTRF, 

(rocRS -^ rocRs^rccRs; ^(0 ^ f(tX T{t)), (3) 

where the transformation matrix T{t) = Q(t)'R(t)'W(t) with 
2(0, ^(0 and W(t) being 3x3 matrices representing the motion 
of the celestial pole in the GCRS, rotation of the Earth around 
the axis associated with the pole and polar motion, correspond- 
ingly, and the derivatives t(t), f(t) are calculated numerically. 

Displacements of station reference points. Station posi- 
tions in the GCRS should be subsequently corrected for vari- 
ous geophysical eff'ects. These include displacements caused by 
the solid Earth tides arising from the direct eff'ect of the external 
tide-generating potential, ocean tidal loading, diurnal and semid- 
iurnal atmospheric pressure loading, and the centrifugal pertur- 
bations caused by the pole tide-driven Earth rotation variations. 
These eff'ects are taken into account following the lERS conven- 
tions 2010 (IPetit& Luzum 2010b . 

Geocentric to Barycentric transformation. Finally, the 
Lorentz transformation needs to be applied to the station GCRS 
position rgc, velocity Tgc and acceleration i^gc to transform it in 
the solar system barycentric reference frame BCRS: 



Tbc 



Tbc 







yUE 



(i+r)c/E 



1 

2? 



(Ve • rgc) • Ve + Re 



(4) 



y2 

2c2 



(1 

+Ve • (1 - ^(Ve • r,e)) 



lgC> 



igc 



(5) 



-^(VE-i^gc)-(VE + 2rge) + AE, 



gc 



(6) 



GMj 



where Lc = 1.48082686741 x 10-^ L^e = Z ^^ is the 



j^i 



TEj 



Newtonian potential of all solar system bodies excluding the 
Earth, evaluated at the geocenter, y is the PPN parameter (equal 
to 1 in General Relativity), Re, Ve and Ae are the position, ve- 
locity and acceleration of the Earth in the BCRS, respectively, 
and c is the speed of light in a vacuum. 
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Time-scale transformations. The time-scale used for mea- 
surement timing at stations is the coordinated universal time 
(UTC), which differs from the physically realised international 
atomic time scale, TAI, by a number of leap-seconds (currently, 
34 seconds). The latter, in turn, differs by 32.184 seconds from 
the terrestrial time (TT), which is the theoretical time-scale for 
clocks at sea-level. The TT is a scaled version of the geocentric 
coordinate time (TCG, the time coordinate of the lAU space- 
time metric GCRS), which eliminates the combined effect on the 
terrestrial clock of the gravitational pote ntial from the Earth an d 
the observatory's diurnal speed in TCG JPetit & Luzum 2Q10l) : 



TCG = TT + Lq- (JDtj - TTo), 



(7) 



where TTq is TT at 1977 January 1.0 TAI, JDtt is TT, both 
expressed as Julian date, and Lq = 6.969290134 x 10"^^. 

The time-scale used in the ephemerides of planetary space- 
craft, as well as that of solar system bodies, is the bary centric 
dynamical time (TDB), a scaled version of the bary centric coor- 
dinate time (TCB) - the time coordinate of the lAU space-time 
metric BCRS. The TDB stays close to TT on the average by 
suppressing a drift in TCB due to the combined effect of the ter- 
restrial observer's orbital speed and the gravitational potential 
from the Sun and pla nets by applying a linear transformation 
(IPetit& Luzum 20101) 



TDB = TCB - Lb- (//)tcb - To)- 86400 -h TDBq 



(8) 



TCB 

Lb = 



expressed as Julian date, Tq 
1.550519768 x 10"^ and TDBq 



where JDjqb is 
2443144.5003725, 
6.55 X 10-5 s. 

The difference between TCG and TCB involves a transfor- 
mation accounting for the orbital speed of the geocenter and the 
gravitational potential from the Sun and planets with the accu- 
racy of 0(c~^) given by 



(9) 



TCB - TCG = c-2( r [ Y + ^E(RE)]d^ + Ve- (R - Re)), 

where Re denotes the BCRS position of the Earth and R is the 
BCRS position of the observing station. 

Computation of the difference TCB - TCG is per- 
formed on the basis of the Fairhead-Bretagnon model 
( [Fairhead & Bretagnon 1990| ) in its full form, which yields on 
accuracy of a few nanoseconds. The full time-scale transfor- 
mation chain is as follows: UTC -^ TAI -^ TT ^ TCG -^ 
TCB -^ TDB. 

VLBI delay model for extragalactic sources. For calcu- 
lating the VLBI signal delay for the extragalactic sources, we 
use the so-called consensus model recommended by the lERS 
(IPetit & Luzum 20101) , which is based on the plane- wave ap- 
proximation when the distance to a source is assumed to be infi- 
nite. In the BCRS, the vacuum delay equation is 



T2-Ti = --• (R2(r2) - Ri(ri)) + ^T., 



(10) 



The total barycentric vacuum delay is then converted into 
the geocentric one by applying the Lorentz transformations and 
is given by 



h-h 



1-h 



K-(VE + r2,gc) 



Ar, 



grav 



Kb 

c 



1- 



(i+r)^E 



2c2 



^E • r2,gc 



VE-b ^^ K-Ve, 



(11) 
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Fig. 2. Geometry of VLBI observations of spacecraft in the 
Barycentric celestial reference frame. 

Here K is the unit vector from the barycentre to the source in 
the absence of gravitational or aberrational bending, b = r2(^l)- 
rl(^l) is the GCRS baseline vector at the time of signal arrival at 
the first station. The general relativistic delay, AJgrav, is given 
for the /th gravitating body by 



Ar. 



grav,/ 



2 — ^m ■ 



/v2/ + -IV' R' 



(12) 



2/ 



with Riy = Riy(ri) = Ri(^i) - R/(^iy) and R2/ = R2/(^i) = 
^2(^1) - "^(K- b) - R/(^i/), where Ri and R2 are the barycentric 
radius-vectors of the first and second receivers, ti is the TCG 
time of arrival of the radio signal at the first VLBI receiver, tu - 
the moment of closest approach of the photon to the gravitating 
body J in TCG. The total gravitational delay is the sum over all 

gravitating bodies including the Earth: AJgrav = Z AT^grav,/ • 

/ 
Delay model for VLBI observations of spacecraft. 

For spacecraft VLBI observations, a delay model should be 
able to account for the "near-field" eff'ects caused by finiteness 
of the source's distance. Indeed, the diff'racted radio wave is con- 
sidered to be in the near field if the distance to the source from 
the aperture R< D^/A (Bom & Wolf 2002 ). where D is the char- 
acteristic size of the aperture and A is the wavelength. For a 
VLBI array with a synthesised aperture of ~ 1000 km observing 
at X-band (A - 3 cm) this yields R < 100 AU. Note that near- 
field delays can diff'er from "classic" ones by microseconds. 

The signal delay in the TDB -frame is defined as the diff'er- 
ence between two light travel times LTi and LT2 from the space- 
craft to the first and the second station of the baseline (Fig. [2]). 
Given the reception time Ti in TDB at the first station, the trans- 
mission time To is calculated by iteratively solving the light-time 
equation for the 'down-leg' 0-T. 



Ti-To = — ^RLTou 

c 



(13) 



where RLTqi is the relativistic term, which accounts for the 
eff'ects of both special and general relativity ( [Moyer 2003 1 ): 



RLToi = 



(l+r)-GMs ^^/?g+7?^+/?gi + 



(l+y)-GMs 

r2 



pS ■ i?S i?S , (l+y)-GMs 



10 






J?B _i_ j?B jpB ' 



(14) 



Here S denotes the Sun and B runs over all solar system 
planets and the Moon, Rq is the spacecraft barycentric positioqj 



Ri 



01 



Ri(ri)-Ro(ro) 



(15) 



^ for the ESA spacecraft, the ephemerides are provided by the ESA 
European Space Operations Centre in Darmstadt, Germany (ESOC) 
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Fig. 3. Kolmogorov spectrum of phase scintillations. 



Rf = Ri(Ti) - R,(Ti), i = 0,Ua = S,B 



(16) 
(17) 



At each iteration step the correction to the Tq estimate is 
given by 



ATn 



Ti-To-Roi/c-RLT^ 



01 



1 -Poi/c 



where 



Poi 



^o(Tol 



(18) 



(19) 



Given the solution Tq of equation ( 13 ), the light-time equa- 
tion for the down-leg 0-2 is solved for reception time T2 in the 
same manner as described above. 

The difference T2 - Ti represents the VLBI delay in the 
bary centric TDB -frame, which is then transformed into the geo- 
centric TT- frame: 



ti-h 



(?f^i4(l-^)i-¥}(-^)""« 



The performance of the geometric delay model described 
here was checked by comparing it with one of the most accurate 
up-to-date analytical models by Sekido & Fukushima (2006), 
which proved to be consistent at the picosecond level. 

Propagation and other additional impacts on the sig- 
nal delay. Random fluctuations of a radio wave phase caused 
by propagation efl'ects are dominated by the interplanetary 
scintillations in the solar wind. The power spectrum of these 
fluctuations is well-represented by the Kolmogorov spectrum 
( ,Kolmogorov 1991j . The scintillation phase is given by 






'/max 

Wo)- f 

fo JO 



f \ slope V'^ 



(21) 



where /)(/o) is the spectral power density, the lower limit of inte- 
gration /o = l/r with T being the length of the phase-referencing 
nodding cycle, the upper limit /max is determined by the system 
noise level and slope is the linear slope of the spectrum in log- 
arithmic scale. For example, with a nodding cycle of r = 180 s, 
/o = l/r ^ 5.5 mHz and the scintillation phase is proportional to 
the marked area in the Fig.|3] Obviously, the shorter the nodding 
cycle, the better the target phase calibration. 

The Earth's ionosphere and troposphere also induce an ad- 
ditional delay to the radio signal propagation. To mitigate the 
efl'ect of the former medium, we use the vertical total electron 
content (vTEC) maps provided by the international GNSS ser- 
vice (IGS) on a daily basis with a two-hour temporal resolution 



on a global grid dFeltens and Schaer 1998b . The IGS uses the 
single thin ionospheric layer model, therefore the vTEC should 
be properly mapped onto the direction of the target source to 
obtain the slant total electron content (TEC): 



TEC 



vTEC 
cosz' ' 



(22) 



where z' is the zenith distance of the target as seen from the 
model's ionospheric layer: 



R 



z = arcsm 



R-^H 



• sm z. 



(23) 



with R being the mean radius of the Earth, H the height of the 
ionospheric layer and z the apparent vacuum zenith distance of 
the target. 

The ionospheric delay for each station can then be computed 
as 



5.308018- r^C 

47r2./2 ' 



(24) 



where / - is the observational frequency (e.g. 8.4 GHz for obser- 
vations in X-band). 

There are two options for mitigating the tropospheric efl'ect 
on the signal propagation delay. They are based either on "ready- 
to-use" empirical models or ray-tracing through the numerical 
weather models (NWM). The former provides the tropospheric 
delays in the direction of the local zenith and a mapping func- 
tion to project the delay onto the directio n to the source. W e use 
the Vienna mapping functions VMFl (B oehm et al. 20061) , rec- 
ommended in the lERS conventions 2010. The precision of the 
empirical models is sufficient for many cases, whereas we em- 
ploy the latter option when an extreme accuracy is required, e.g. 
if the target sources is very low above the horizon and/or the an- 
gular distance between the spacecraft and the phase-referencing 
source is large. The latter was the case in the VEX tracking ex- 
periment EM081c, which we describe in the next section. We use 
the ray- tracing algorithm from Hobiger et al. C2008J with some 
modifications adopted from Duev et al. (12011b . The meteorolog- 
ical data from the European centre for medium-range weather 
forecasts (ECMWF) are used. 

Additionally, contributions to the delay due to the an- 
tenna axis off'sets and smaller size eff'ects due to ther- 
mal ([Nothnagel 2009|) and gravitational dSarti et al. 20091 
'Clark & Thomsen 1 988) deformations of telescopes are also 
taken into account in our delay model. 

2.6. Estimation of corrections to the S/C state vector 

The main task of VLBI tracking of spacecraft is estimating the 
state vector as a function of time. This task has some similari- 
ties with traditional astrometric applications of the VLBI tech- 
nique. However, because a spacecraft is in the near-field of the 
VLBI-synthesised aperture, several modifications of traditional 
approaches and formalisms of far-field VLBI are necessary. 

Analogues of t/v-coverage for near-field VLBI. The uv- 
projections of baselines used in the standard far-field VLBI are 
basically the degenerated (in the limit of an infinite distance to 
the source) partial derivatives of the predicted geocentric delays 



by the source position lateral to the line of sight: 



■ dipTij' c. 



^ deTi, 



(25) 



where c is the speed of light in a vacuum and ip and 6 are the 
geocentric spherical angular coordinates of the source. 
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For near-field VLBI observations at the epoch t the uv- 
projections of the baselines i-j (i=l,N-l,j = i-\-l,N, where 
N is the number of participating stations), because of the finite- 
ness of the distance to the source, should be replaced with the 
corresponding elements of the Jacobian /: 



(26) 



where t/ is the geocentric delay for the station / at t and 
the partials in the Jacobian are computed numerically. In other 
words, the traditional wv-formalism is merely a special case (a 
target source at an infinitely large distance) of a more general 
formalism, the Jacobians. 

Spacecraft imaging and state vector estimation. 

The ('dirty') image of the spacecraft at the epoch t is recov- 
ered by performing the Fourier transform of the sampled "visi- 
bility" function (see [Thompson et al. 200 1| ) 



/ dm 


dT]-> \ 

de 




/ d(Ti-T2) 

dip 


d(Ti-T2) \ 

de 


dTiN 

dip 


dTiN 

de 


= 


d(Ti -tn) 
dip 


d(Ti -Tn) 

de 


dTN-l,N 


dTN-l^N 




d(TN-l-TN) 


d(TN-l-TN) 


^ dip 


de ) 




^ dip 


de ) 



ItiUm) 



= f%{S(u, 



v)' V(u, v)' e 



-2m-{ul+vm) 



) dw dv , 



(27) 



with the sampling function being the sum over all baselines 
of two-dimensional delta functions: 

S (u, ^) [ = ^ ^ ^(^ - ^ij^ ^ - ^ij) [ (28) 

and the normalised weighted 'visibility' function given by 

-2m-(pii 



V(uij, Vij) = w/ 



(29) 



where Wfj is the weight coefficient for the baseline i-j at time t. 

For a given time range t = [tsJe], the resulting image is 
reconstructed as a sum over all epochs within the time range: 



1(1, m)= ^/,(/,m) 



A^-l A^ 



" Z Z Z ^^J' 5l(e^-2^-^"'^-^^^'^'^-'^'^-^) 
The measurement equation for the epoch t is 



(30) 



(31) 



where Acp is a [A^b x 1] vector of diff'erential phases on all base- 
lines (A^b = N' (N- 1)/2 is the number of baselines) and Aa is the 
vector of corrections to the a priori geocentric angular position 
of the spacecraft: 



A0 = 



^12 



^lA^ 



Aa : 



_(A<p 
AO 



(32) 



V (PN-l,N J 



Finally, the astrometric solution of equation ( [3T] ) in the least- 
squares sense for each epoch t is given by 



a2| =((/^-/)-^-/^-a0)|. 




(33) 
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X, arcsec (from west to east) 

Fig. 4. VEX orbit around Venus on 201 1 March 28, as seen from 
the geocentre. Red boxes represent the position of VEX at the 
epochs when imaging was performed. 

The main purpose for the introduction of the Jacobian for- 
malism is to use it for the estimation of the corrections to the 
a priori position of a near-field VLBI target according to equa- 
tions ([3T])-([33]). Imaging, although it can be used for the same 
purpose yielding the same result, is merely a "by-product" of 
such a formalism, and is used in our paper mostly for the pur- 
poses of illustration. In addition to this, the closer the target is 
to an observer, the more significant the diff'erence between the 
Jacobians and "normal" wv-projections becomes. 

3. VEX VLBI observations and results 

ESA's Venus Express Spacecraft has been observed with the 
European VLBI Network (EVN) within the framework of the 
PRIDE initiative under the EVN observing project EM08L 

From four requested observing epochs we used only three. 
The first two observing runs, on 2010 August 23 (2 hours long, 
3 stations) and on 2010 September 20 (3 hours long, 2 sta- 
tions) were mainly devoted to tuning-up and verification of the 
technique, namely developing and testing the S/C observation- 
related scheduling, data acquisition, transfer and processing 
methods with the near-field delay model. 

Here we report the results obtained with the third obser- 
vational run, which has been the most intensive experiment 
so far in terms of the number of participating telescopes: the 
EVN stations Onsala (On, Sweden), Wettzell (Wz, Germany), 
Medicina (Mc, Italy), Matera (Ma, Italy), Yebes (Ys, Spain), 
Metsahovi (Mh, Finland), Svetloe (Sv, Russia), Zelenchukskaya 
(Zc, Russia), Hartebeesthoek (Hh, South Africa), and the NRAO 
Very Long BaseHne Array (VLBA) station St. Croix (Sc, USA). 

The observations were conducted on 2011 March 28, 
from 08:45 to 11:30 UTC. The source J2225-0457 (also 
known as 3C446) was used as a fringe finder and for the 
clock search. The source J221 1-1328 with J2000.0 coordinates 
of 22^1^24^.0994590, -13°28'09.723950'' (Petrov2011) was 
used as a phase-reference calibrator. The estimated coordinate 
error for this source is 0. 1 1 mas in Right Ascension and 0. 12 mas 
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Fig. 5. Near-field Jacobians for the VEX pointings in kilometers, 
201 1 March 28, the EVN project EM081. The longest east- west 
baseline - Ys-Zc, the longest north-south baseline - Mh-Ys 



in Declination JPetrov 201 lb . All involved EVN radio telescopes 
were configured to observe the spacecraft signal in a phase- 
referencing mode, switching between the S/C and the phase ref- 
erence source with a nodding cycle of 4 minutes, including a 
20 second gap for re-pointing and antenna calibration. At On, 
Mc, Ma, Mh, Hh, we recorded 79 scans per station, including 
2 for the fringe finder, 39 for the calibrator source and 38 for 
VEX. The stations Ys and Sc joined in starting from the 15^^ 
scan owing to technical problems and limited visibility, respec- 
tively. The eastern-most stations Sv and Zc stopped observ- 
ing according to schedule owing to the target sources' setting. 
During the observing run, the VEX spacecraft had mean co- 
ordinates 22^U'^l3\-ir4V2T\ At the time of the observa- 
tions, Venus was at a solar elongation of 36 degrees and at a 
distance from Earth of 1.24 AU. Fig. |4] shows the orbital posi- 
tion of VEX around Venus during the observations as seen from 
the geocentre. Fig. [5] shows the pseudo-wv-coverage (Jacobians) 
of the EM081C VLBI array. 

It should be pointed out that the observational setup for the 
experiment was not ideal - target declinations ranging from 
-11° to -13° are unfavourable for the EVN array. In addition 
to this, the angular separation of ~ 2.5° between Venus and the 
calibrator was at the edge of the isoplanatic patclj^ at low el- 
evations caused by a low declination of the target, which lead 
to a relatively large unmodelled tropospheric/ionospheric phase 
diff'erence between the calibrator and VEX. 

The data were acquired using the standard VLBI data ac- 
quisition system based on the Mark5 A/B recording equipment. 
We used a base frequency of 841 1.99 MHz, four 16-MHz-wide 
channels with 2-bit Nyquist sampling, resulting on aggregate 
data rate of 256 Mbps per station. The data files were either 
transferred over the high-capacity network connection after the 
experiment to Metsahovi and JIVE or recorded on disks and then 
shipped to the processing centre at JIVE. 



^ Isoplanatic patch is an angular region on the sky within which the 
measured phases of a calibrator can be coherently applied to the source 
phases. Its size is determined by the atmospheric turbulence. 
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Fig. 6. Frequency Doppler shift over the path of the ES A track- 
ing station Cebreros - VEX - Metsahovi, Medicina, Yebes (top 
left), Cebreros - VEX - Metsahovi, Medicina, Yebes - geocentre 
(top right). The diff'erential frequencies in the geocentre with a 
standard deviation of- 5 mHz at 2.5 second sampling are shown 
at the bottom. 2010 August 23, EVN project EM081. 




Fig. 7. Interferometric fringes of the VEX signal in the 
spectral domain, 1024 frequency points, baseline Onsala- 
Zelenchukskaya, 2011 March 28, EVN project EM081, using 
the near-field delay model and two-way Doppler phase correc- 
tion to avoid frequency smearing when correlating at the SFXC 
at JIVE. Top: amplitude in relative units; bottom: phase in radi- 
ans. The horizontal axis - frequency in video band (MHz) rela- 
tive to 8.41 1990 GHz. Most of the power of the spacecraft signal 
is concentrated within a ~1 MHz- wide band around the carrier 
line, although the S/N for the carrier is considerably higher than 
for the rest of the band. 



The narrow-band data processing and carrier fre- 
quency/phase extraction were conducted with the SWSpec 
and SCtracker software at Metsahovi. The Doppler frequency 
detections from all stations reduced to the geocentre using 
the near-field delay model (as /gc(0 = ftcO - t)'(^ - t)) 
diff'er on average only by several mHz or even less, proving 
the consistency between the actual and model delays. As an 
example, the topocentric frequency detections measured for the 
stations Mh, Mc and Ys, detections reduced to the geocentre and 
diff'erential geocentric frequencies for the experiment conducted 
on 2010 August 23 are shown in Fig. [6] 

The broadband correlation was conducted at JIVE using the 
SFXC. The fringe finder data were used to perform a clock 
search. The VEX cross-correlation spectrum is smeared due to 
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Fig. 8. (a, b) - Fringe finder (3C446, 2.5 Jy source at an angular distance of 15° north from the target) and a calibrator (J2211- 
1328, 0.17 Jy source at an angular distance of 2.5° south from the target) self-calibrated images using the EVN array excluding 
Hh, (c, d) - self-calibrated images of the same sources obtained with the whole array. Horizontal axis - deviation from the nominal 
source's Right Ascension in mas, vertical axis - deviation from the nominal source's Declination in mas. Dashed white ellipses at 
the bottom-right represent the synthesised beam size and shape. 




Fig. 9. Radio images of VEX obtained with the EVN array. The horizontal axis represents the deviation from the nominal source's 
Right Ascension in mas, vertical axis - deviation from the nominal source's Declination in mas. (a) 09h05m TDB, stations On, Wz, 
Ma, Mc, Mh, Sv, Zc, (b) 09h30m TDB, Ys joined in, (c) 09h55m TDB, Sv left, (d) 10h20m TDB, Mh left, (e) 10h45m TDB, Zc 
left. No cleaning applied. 



the intrinsic change of the frequency emitted by the S/C, as it re- 
transmits the signal in the two-way Doppler regime. This change 
may reach several kHz over a time span of one hour, which is 
shown in Fig.|6] To avoid this smearing when correlated with the 
SFXC, a two-way Doppler correction (see formula (T) and the 
paragraph below it) was applied to the S/C data. Fig.JTjrepresents 
the amplitude and phase of the cross-correlation function of the 
VEX signal (amplitude and phase) on the baseline On-Mh. Most 
of the spectral power is concentrated within a ~ 1 MHz region 
around the carrier line, allowing us to resolve the In ambiguity 
between the group delay detected by the broad-band correlation 
with the SFXC and phase delay detected using the SCtracker 
software. Fig.[8](a,b) shows the self-calibrated and cleaned im- 
ages of the fringe finder and the calibrator obtained with the 
Difmap ( [Shepherd 1995] ) and MPS dGreisen 2003 ) software. 

Unfortunately, an anomalous clock rate off'set at Hh did not 
allow us to obtain its absolute phases. Although good for self- 
calibrated imaging (see Fig. [S]), the Hh data were not used for 
absolute positioning during the EMOSlc observational run on 
2011 March 28. The Hh LO converter was fixed several days 
later, based on our analysis of the LO behaviour made with ultra- 
high spectral resolution SCTrack software package. 

The correlation of the calibrator source for each baseline pro- 
vides the residual delays. For coherent phase-referencing, the 
residual delays must be lower than the radio wave period at X- 
band, which is about 0.120 ns. The measured residual delays and 
phases were applied to the detected residual phases of the target. 
With these we performed imaging of the VEX spacecraft and 
derived the deviations of the lateral coordinates of the spacecraft 
with respect to a priori ones. Fig. [9] shows the imaging results ob- 
tained for VEX with the EVN array without any cleaning. The 



averaging time for each image is 20 minutes. Fig. ^Q| shows the 
estimated deviations of the a priori lateral position of the VEX 
spacecraft. The integration time for each point is 30 seconds. The 
coordinate off'sets were estimated using equation 33 However, 
for a consistency check, we made S/C images for each short in- 
terval of 30 seconds with a "super"-resolution and measured the 
peak position on each of them. The results yielded a consistency 
of the two approaches at the level of 25 microarcseconds. The 
orbit accuracy estimate is at a 3 sigma level of 200-300 m across 
the track and 500-600 m along the track. 

The north/south VEX off'sets of ~ 1 mas from the expected 
orbital position seen in Fig. [T0| although nearly within the er- 



ror budget, may be caused by systematic phase errors between 
the target and calibrator due to their low elevation. More exper- 
iments with diff'erent relative positions of the S/C and calibra- 
tor(s) are needed to prove or disprove this claim. 

4. Discussion and conclusions 

The EVN experiment EM08 1 enabled us to show that astromet- 
ric spacecraft positioning with a very high accuracy is achiev- 
able using the PRIDE approach. A relatively large angular dis- 
tance between VEX and the calibrator source (2.5°), their low 
elevation and a limited accuracy of the troposphere/ionosphere 
models did not allow us to include the data from St. Croix (Sc) 
in the astrometric solution without a baseline fringe phase In- 
bifurcation. Fig.[8](c, d) shows the self-calibrated images of the 
fringe finder and the calibrator with Sc and Hh added. As seen 
from this figure, more than a three-fold resolution improvement 
(at the level of up to 20 mas, 0.1 nrad) is possible. We adopt this 
value as the aim for our future experiments. To achieve this in- 



D.A. Duev et al.: S/C VLBI & Doppler tracking: algorithms and implementation 



F 


? 


^ 




<n 








Q) 






1 


it 




O 




> 




o3 
X 






,; w: 



\v 






r%\ 



• ~» 



"*< 









1^% 



' • * 




11.5 



-2 -1 1 
X offset, km 

Fig. 10. (top) Measured coordinate offsets with respect to the a priori trajectory, supplied by ESOC. The X-offset is shown in blue, 
Y-ofFset - in red. X and Y are approximately RA and Dec. The measurement sampling - 30 s, 1 km corresponds to ^ 1.1 mas at the 
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crease in resolution, we will have to resolve the In ambiguities 
for the longer baselines. For this purpose, closer (in-beam in the 
ideal case), preferably multiple calibrators should be used in the 
experimental setup for phase referencing. In addition to this, bet- 
ter troposphere and ionosphere models would be very important, 
especially for the low elevations of the sources. 

The results of the current study in terms of achieved accuracy 
are compatible with those of the other recent well-established 
VLBI spacecraft tracking projects, such as the VLBA tracking 
of the Cassini spacecraft (J Jones et al. 20111) and NASA's Mars 
Exploration Rove r B spacecraft during its final cruise phase 
( |Lanyi et al. 2007| ). 

The PRIDE approach described in this paper proves its appli- 
cability to virtually any deep-space mission almost anywhere in 
the solar system. We note that the approach is not exclusively 
oriented on planetary science missions, but other deep-space 
and near-Earth targets. This has been d emonstrated by VLBI 
observations of GL QNASS satellites ("Tomat ore et al. 20101 
ITornatore et al. 20lT1) aiming to establish a better link between 
the terrestrial and celectial reference frames. 
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